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Abstract. We present evolutionary models for cool brown dwarfs and extra-solar giant planets. The models re- 
produce the main trends of observed methane dwarfs in near-IR color-magnitude diagrams. We also present 
evolutionary models for irradiated planets, coupling for the first time irradiated atmosphere profiles and inner 
structures. We focus on HD 209458-like systems and show that irradiation effects can substantially affect the ra- 
dius of sub-jovian mass giant planets. Irradiation effects, however, cannot alone explain the large observed radius 
of HD 209458b. Adopting assumptions which optimise irradiation effects and taking into account the extension 
of the outer atmospheric layers, we still find ~ 20% discrepancy between observed and theoretical radii. An extra 
source of energy seems to be required to explain the observed value of the first transit planet. 
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1. Introduction 

The past decade was marked by two major discoveries 
in the field of stellar and planetary physics: the detec- 
tions of the first unambiguous brown dwarf (BD) GL 229B 
(Oppenheimer et al . 1995) and the first extrasolar giant 
planet (EGP) 51 Peg b (Mayor & Queloz 1995). The near- 
IR spectrum of GL 229B was found to be dominated by 
strong methane absorption bands, looking more similar 
to Jupiter than to late type-stars. On the other hand, the 
surprisingly small orbital separation between 51 Peg b and 
its parent star suggests that the planet should be affected 
by irradiation and that, given the expected large surface 
temperature, its atmospheric properties should resemble 
more the ones of relatively hot brown dwarfs than the ones 
of jovian planets. 

Since then, about thirty methane dwarfs (or the so- 
called T-dwarfs) have been identified, due mainly to 
the near-IR surveys 2MASS (Burgasser et al. 1999), 
SDSS (Strauss et al. 1999) and the VLT (Cuby et al. 
1999). The radial velocity technique has now revealed 
more than 100 EGPs in orbit around nearby stars (see 
Hubbard, Burrows, Lunine, 2002 for a review and ref- 
erences therein), with a large fraction (~ 10%-20%) be- 
ing extremely close (less than 0.06 AU) to their parent 
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star. The mass of substellar companions detected by ra- 
dial velocimctry extends well above the deuterium burning 
minimum mass O.O12M0 = 12Mj (Saumon et al., 1996; 
Chabrier et al., 2000a). This mass is often used as the 
boundary between planets and brown dwarfs, more for 
semantical than physical reasons. On the other hand iso- 
lated objects with planetary masses are now discovered in 
young stellar clusters, down to a few Mj, as recently re- 
ported by Zapatero et al. (2002) in a Orionis. These obser- 
vations suggest that there is an overlap between the mass 
range of the least massive brown dwarfs and of the most 
massive giant planets. In principle, different formation 
processes should distinguish planets from brown dwarfs. 
However, such a distinction is difficult to characterize in 
terms of atmospheric, structural and cooling properties 
since both types of objects have convective interiors with 
essentially a metallic H/He mixture. The signature of a 
central rock+ice core, like in solar giant planets, would be 
the clear identification of a planet. The presence of a core 
can affect the radius of a planet, yielding a smaller plan- 
etary radius than in the absence of a core. For 1 Mj, the 
effect is about 5% on the radius for a core mass < 0.06 
Mj (see Saumon et al. 1996). The presence of this core 
can be inferred from the accurate characterization of the 
gravitational moments of the object, and such an observa- 
tion is currently not feasible for EGPs. In addition, both 
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giant planets and brown dwarfs have atmospheres domi- 
nated by molecular absorption and the effects of cloud for- 
mation. Although frustrating from an observational point 
of view, these similarities imply that the general cooling 
theory for BDs, involving detailed models of the atmo- 
sphere and inner structures, can be applied to EGPs. In 
terms of cooling properties, this general theory can even 
be applied to Jupiter, as emphasized by Hubbard et al. 
(2002). Additional observational constraints, as provided 
by spacecraft encounters or by direct probes (e.g Galileo), 
have lead to refinements of the models (heavy clement 
core, non-standard chemical composition). As mentioned 
above, such constraints arc, unfortunately, far from being 
accessible for EGPs. 

Much effort has been devoted to the modeling of sub- 
stellar objects during the past decade, improving our un- 
derstanding of cool atmospheres (see Allard et al. 1997 for 
a review), of the role of dust (Tsuji et al. 1999; Burrows 
et al. 2000; Ackerman & Marley 2001, Allard et al. 2001; 
Marley et al. 2002), of irradiation (Saumon et al. 1996; 
Seager & Sasselov 1998; Sudarsky et al. 2000; Barman et 
al. 2001), and of their inner structure and evolutionary 
properties (Burrows et al. 1997 ; Chabrier et al. 2000b; 
see Chabrier & Baraffe 2000 for a review) . One remaining 
major challenge in the theory is the description of dynami- 
cal processes of grain formation and diffusion necessary to 
understand the transition between L-dwarfs and T-dwarfs, 
which is expected to take place at T c g ~ 1300K - 1700K. 
The former objects are better reproduced by dusty atmo- 
sphere models, whereas the later are better reproduced 
by dust-free (or partly dusty) models. The recent observa- 
tions of L/T dwarfs at the transition clearly indicate that 
complex processes take place in the atmosphere of these 
objects (see e.g Burgasser et al. 2002). Another important 
challenge is the modeling of irradiation effects, which are 
expected to affect the spectra of close-in EGPs, and may 
also affect their inner structure and cooling properties. 
The recent discovery of the planet HD209458b transit- 
ing its parent star (Charbonneau et al. 2000) provides a 
unique test to explore such effects, since its mass and ra- 
dius can be determined with high accuracy from the mod- 
eling of the transit lightcurve. According to the most re- 
cent determination (Cody & Sasselov 2002), the mass and 
radius of the planet are estimated to be m — 0.69±0.02Mj 
and R = lA2toi° 3 Rj. 

Evolutionary models including crude estimates of the 
effects of irradiation on planet atmospheres suggest that 
extrinsic heating is sufficient to maintain a larger plan- 
etary radius compared to an isolated planet. It has thus 
been argued that irradiation could explain the large radius 
of HD209458b (Guillot et al. 1996; Burrows et al. 2000). 
More recently, Guillot & Showman (2002) questioned such 
results and argue that the radius of HD209458b can only 
be reproduced if the deep atmosphere is much hotter than 
what can be expected from irradiation effects. However, 
none of these calculations includes a consistent treatment 
between the irradiated atmospheric structure and the inte- 
rior structure of the planet. Such a consistent treatment 



is mandatory to get reliable results since the deep interior 
entropy profile, which determines the heat content of the 
planet to be radiated away while it cools, is affected by the 
modification of the atmospheric temperature profile due 
to the incoming external heat flux. The main goal of the 
present paper is to present the first such consistent calcula- 
tions. As mentioned above, in the case of non-irradiation, 
these calculations apply to the evolution of cool (dust- 
free like) brown dwarfs, i.e. T-dwarfs, and extrasolar giant 
planets far enough from the parent star for the irradiation 
effects on the thermal structure of the planet to be negligi- 
ble. This is the case of the solar giant planets, the cooling 
of which is simply characterized by the cooling properties 
of the " isolated" planet plus the heating contribution from 
the Sun AitctR^Tq, where R p is the radius of the planet 
and Tq represents the equivalent black body temperature 
of the converted solar radiation (Hubbard, 1977; Guillot et 
al., 1995). In section 2 we briefly present the input physics 
of non-irradiated models, describing methane dwarfs and 
isolated EGPs. Apart from the impinging stellar flux, the 
same input physics are used to analyse the effects of ir- 
radiation. The effects on the radius and cooling proper- 
ties of giant planets are described in section 3 and results 
are compared to the observed properties of HD209458b. 
Discussion follows in §4. 

2. Non-irradiated models 

2.1. Model description 

The main input physics involved in the present models are 
the same as described in our previous works (Chabrier & 
Baraffe 1997; Baraffe et al. 1998; Chabrier et al. 2000b). 
The models are based on the coupling between interior and 
non-grey atmosphere structures. The treatment of dust 
in the atmosphere is described in detail in Allard et al. 
(2001), with two limiting cases of dust treatment. The 
first case, referred to as "DUSTY" , takes into account the 
formation of dust in the equation of state, and its scat- 
tering and absorption in the radiative transfer equation. 
Such models assume that dust species remain where they 
form, according to the chemical equilibrium conditions. 
The second case, referred to as "COND", neglects dust 
opacity in the radiative transfer equation. In a previous 
paper (Chabrier et al. 2000b), we presented the evolution- 
ary models based on DUSTY atmosphere models, aimed 
at describing the evolution and the photometric and spec- 
troscopic properties of early L-dwarfs. The present paper 
is devoted to evolutionary models based on the COND 
approach, which are more appropriate to objects with ef- 
fective temperatures T e s < 1300K, such as methane dwarfs 
or EGPs at large orbital separation. These models apply 
when all grains have gravitationally settled below the pho- 
tosphere. A forthcoming paper will be devoted to models 
taking into account characteristic diffusion timescales of 
different processes affecting the dust stratification (e.g co- 
agulation, gravitational settling, convection). These mod- 
els aim at describing in particular the transition objects 
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between late L-dwarfs and early T-dwarfs (see Allard 
2002). 



2.2. Evolution of cool brown dwarfs 

A preliminary version of the COND models was presented 
in Chabrier et al. (2000b) down to 0.01 M©. In the present 
paper, we extend our calculations to T e g = 100 K and m = 
0.5 Mj. The evolution of L and T c g as a function of time 
for different masses is displayed in Fig. ^ The properties 
of the COND models for different ages are given in Tables 
1-5. As already stressed in Chabrier et al. (2000b), the 
treatment of dust in the atmospheric models barely af- 
fects the evolutionary tracks in terms of L and T e g as a 
function of time for a given mass (see Fig. 2 of Chabrier et 
al. 2000b). Consequently, although the COND models are 
more appropriate to describe the spectral and photomet- 
rical properties of substcllar objects with T ff < 1300K, 
they provide a good description of their cooling properties 
even at higher T e g. In other words, it is not necessary to 
perform evolutionary calculations with the DUSTY mod- 
els above T e g ~ 1300K and switch to the COND models 
below. An object characteristic of the present models was 
recently discovered by Zapatero et al. (2002) in a Orionis: 
S Ori 70. From a comparison of its observed spectrum with 
COND synthetic spectra (Allard et al. 2001), Zapatero et 
al. (2002) estimate an effective temperature T c s ~ 700- 
1000K. If the membership of S Ori 70 to a Orionis is con- 
firmed, implying an age < 10 Myr, its mass mass should 
be <5Mj (see Fig.[TJ). 

Several methane dwarfs have been discovered in the 
solar neighbourhood, implying older ages and thus larger 
masses than the extreme case of S Ori 70. At an age of 10 8 
yr, only objects with masses below the deuterium burn- 
ing minimum mass (m < O.O12Af ) have T ff < 1300if , 
whereas at 5 Gyr, it is the case for all substellar ob- 
jects with m<0.06M© (see Fig. ^ and Tables 1 and 
4) . Photometric observations and parallax determinations 
of several L- and T-dwarfs (Els et al. 2001; Leggett 
et al. 2002a; Dahn et al. 2002) now allow a compari- 
son with models in observational color-magnitude dia- 
grams (CMD), providing stringent constraints on theo- 
retical models (see Figs. EE}- 

As already noticed in Allard et al. (1996) for GL 
229B, models free of atmospheric dust clouds better re- 
produce the near-IR photometric and spectral properties 
of methane dwarfs. This is illustrated in (J — K) and 
(K — L') colors in Fig. [2] and Fig. |3| respectively, where 
the COND models reproduce the main observed trends. 
In Fig. EI we note the two transition objects, intermediate 
between L- and T- dwarfs, with (J — iC) ~ 1 (GL 86B: 
Els et al. 1999; SDSS 1254-01: Leggett et al. 2002a, Dahn 
et al. 2002) and the faintest L-dwarfs (Dahn et al. 2002), 
which arc not described by either the DUSTY or COND 
limiting cases, and require a detailed treatment of dust 
diffusion in the atmosphere, as mentioned in §2.1. 
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Fig. 1. Evolution of the luminosity L and effective temper- 
ature T e ff as a function of time (in yr) for different masses. 
Substellar objects are indicated by solid lines and stars 
by short-dashed lines. The horizontal long-dashed line in- 
dicates the limit (T cff ~1300K) below which the COND 
models are appropriate for the photometric and spectro- 
scopic description of T-dwarfs and EGPs (see text). 



The predictions of the COND models provide a gen- 
eral good agreement with observed near-IR photometry 
and spectra at wavelength > 1/xm (Leggett et al. 2002b; 
Zapatero et al. 2002). The models show however short- 
comings at shorter wavelength, with a flux excess around 
0.8-0.9 /xm, characteristic of the I-bandpass. This problem 
is illustrated in Fig. H in a (I — J) - Mj CMD, where the 
COND models predict significantly bluer (/ — J) colors 
than observations. As mentioned in Allard et al. (2001), 
uncertainties in the current treatment of the far wings of 
the absorption lines of alkali elements (Na, K) at such 
pressures may be responsible for this discrepancy. No the- 
ory, however, exists to date for an accurate description 
of broadening of atomic lines by collisions with H2 and 
He. Attempts to improve current treatments are under 
progress (Burrows & Volobuyev 2002). 

The correct trend of colors and spectral properties pre- 
dicted by the present models at wavelength > 1/im, where 
most of the flux is emitted for the concerned range of T e g, 
comfort us however with their reliability to describe ex- 
tremely cool objects. 
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(J-K) 

Fig. 2. Color - Magnitude diagram (J — K) - Mk- 
Observations are taken from Leggett (1992) (mostly for 
M-dwarfs) and Dahn et al. (2002). Also shown: LHS 102B 
(Goldman et al. 1999), GL86 (Els et al. 2001). M-dwarfs 
are shown by dots, L-dwarfs by filled squares and T-dwarfs 
by triangles. DUSTY isochrones (Chabrier et al. 2000b) 
are displayed in the upper right part of the figure, for 
different ages, as indicated. The COND isochrones are 
displayed in the left part of the figure. Some masses (in 
M Q ) and their corresponding T c g are indicated on the 1 
Gyr isochrones by open squares (COND) and open circles 
(DUSTY). The names of two L/T transition objects and 
of the faintest T-dwarf known with parallax are indicated. 

3. Irradiated models 

3.1. Effect on atmosphere structure 

As mentioned in the introduction, a non-negligible frac- 
tion of EGPs orbit close to their parent star and their 
thermal and mechanical structure is affected by irradia- 
tion effects. Therefore, a general theory of cool substel- 
lar objects must take these effects into account. Recently 
Barman et al. (2001) have modeled irradiated atmospheres 
by including the impinging radiation field in the solution 
of the radiative transfer equation. As shown by these au- 
thors, for a given intrinsic luminosity, non-irradiated plan- 
ets have very different temperature structures than irra- 
diated planets. Thus, substituting non-irradiated atmo- 
spheric structures with T e g — T eq (see definition below, 
Eq.(8)) for irradiated structures, as done up to now in the 
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Fig. 3. Color - Magnitude diagram (K — L') - Mk. K is in 
the MKO-NIR system. Observations are from Leggett et 
al. (2002a). Symbols are the same as in Fig. El A DUSTY 
isochrone of 1 Gyr (Chabrier et al. 2000b) is indicated by 
the long-dashed line. The COND isochrones are displayed 
for 0.1 Gyr (dash), 1 Gyr (solid) and 10 Gyr (dash-dot). 
Some masses (in M Q ) and their corresponding T cff are 
indicated on the 1 Gyr COND isochrone by open squares. 



literature, yields incorrect inner boundary conditions for 
evolutionary calculations (see e.g. figure 13 of Barman et 
al. (2001)). Given the present lack of an accurate treat- 
ment of atmospheric dust diffusion, the calculations were 
performed only for the DUSTY and COND limit cases, re- 
spectively. The results emphasize the strong dependence of 
the emergent spectrum and atmospheric structure on the 
presence or absence of dust. In the absence of dust, the 
impinging flux can penetrate in deeper layers of the planet 
atmosphere, affecting more drastically the inner structure 
of the planet than in the dusty case. 

Except for a possible detection of sodium absorption 
in the atmosphere of HD 209458b (Charbonneau et al. 
2002), no constraints on the atmospheric composition of 
EGPs are available at the present time. The only strong 
observational constraint available for irradiated models is 
the transit planet HD 209458b. The determination of its 
mass and radius provides a stringent test to irradiated 
atmosphere calculations and to the resulting structure and 
evolution. We thus apply our calculations of irradiated 
EGPs to HD 209458-like systems. 
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Observations are from Leggett (1992) and Dahn 
(2002). Symbols and curves are the same as in Fig. 



Mj. 
et al. 

m 



We have computed a grid of irradiated atmosphere 
models based on the COND input physics described in 
§2, as in Barman et al. (2001). Although more appropriate 
for EGPs with T cS < 1300 K, the COND models maximise 
the effect of irradiation on the inner atmosphere structure 
and thus on the evolution of EPGs (Barman et al. 2001). 
The grid covers a wide range of T c g from 40K to 100K, in 
steps of 20K, and from 100K to 2800K, in steps of 100K. 
It covers a range of surface gravities from log g = 2.5 to 
log g — 4.5, in steps of 0.5 dex. We adopt the character- 
istics of HD 209458, assuming for the primary an effec- 
tive temperature T e ff* = 6000 K, a radius R± = 1.18 Rq 
(Mazeh et al. 2000; Cody & Sasselov 2002) and an or- 
bital separation a — 0.046 AU (Charbonneau et al. 2000). 
As in Barman et al. (2001), we make the simplifying as- 
sumptions that the impinging radiation field is isotropic 
and the incident flux F inc is redistributed only over the 
dayside, i.e. 



the transfer equation, is the intensity integrated over both 
in-coming and out-going directions fJb{^i = cos 6, where 9 is 
the angle of incidence) . Assuming there is no extra source 
or sink of energy (e.g no horizontal energy transfer), en- 
ergy conservation implies that all the incident energy com- 
ing in must go out. Therefore, in the case of irradiation, 
the in-coming flux from the parent star cancels out the 
extra out-going, absorbed and reradiated flux due to the 
heating of the upper layers of the planet atmosphere (see 
Fig. EJ. The in-coming flux at the surface is F; n = - F{ nc 
and the out-going flux at the surface is F out = F inc + 
crT 4 ff , where crT 4 ff defines the intrinsic, unperturbed flux 
cT 4 ff of the initial, non-irradiated atmosphere structure. 
Energy conservation thus implies: 



Fr>p.t — Foi 



F in = <JT, 



eir 



(2) 



The non-irradiation case (Fj nc = 0) corresponds to the 
usual condition F not = F out = crT 4 ff . 

Our atmosphere models, irradiated or not, are thus 
characterized by the parameters T c g and g. Of course, 
the same net flux F nct corresponds to two different at- 
mospheric structures, in the non-irradiated and irradiated 
case, because of the extra energy source Fj n ^ in the 
latter case (see Fig. 0. Given the above definitions, the 
net flux characterizes the intrinsic luminosity, i.e. the rate 
of energy released by the planet as it contracts and cools 
down: 



Li 



47ri?p(rT c 4 ff = 



—T——dm. 
dt 



(3) 



This quantity determines the cooling properties of the 
planet for a given set of outer boundary conditions pro- 
vided by the atmospheric profile (see §3.2 below). We 
stress that, in the case of irradiation, T e g does not charac- 
terize the total flux emitted by the planet, which is given 
by: 



F ,,t = OT, 



err 



F — rrT 4 

-fine. — O-L^, 



I a 



(4) 



^ t 4 nerm = <fl% + (1 - A)F inc , 



^ine = |(^) 2 F„ (1) 

2 a 

where F* is the total flux from the primary (see discussion 
in §4). 

Before proceeding any further, we briefly re-specify 
definitions of fluxes (see e.g. Brett & Smith 1993), since 
use of various terminologies leads to confusion. In all cases, 
the integrated net flux F ne t, obtained from the solution of F out = crT^ heim 



Note that F ont is the important quantity for observers, 
since it characterises the total radiation of the planet, 
including both thermal and reflected parts of the flux. 
However, we do not focus on this quantity, since a forth- 
coming paper will be devoted to spectral properties of 
irradiated planets (Barman et al. 2003, in preparation). 

For the sake of comparison with non-irradiated atmo- 
sphere profiles (see Barman et al. 2001), we also define 
the quantity T t herm which characterises the thermal flux 
crr t 4 lcrm emitted by the irradiated fraction of the planet (in 
the present case, the day side only). This quantity reads: 



(5) 



where A is the Bond albedo. According to the definitions 
above: 



F re fl, 



(6) 
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where 



F refl = AF ir 



(7) 



is the reflected part of the incident flux. 

Within the conditions of the present calculations 
(Eq(l), a = 0.046, T off * = 6000 K, i?* = 1.18 R Q ), 
our Bond albedo is close to 0.1 for the coolest models 
(T G ff ~ 100K). A final quantity, often used in the litera- 
ture, is the equilibrium temperature, T eq , which character- 
izes the planet's luminosity after having exhausted all its 
internal heat content (see e.g Guillot et al. 1996; Saumon 
et al. 1996): 



4 _ (l-A) p 



2a 



(8) 



T t 4 herm when T c ff -> 0. 



Note that given our definition of F- Inci T c 4 q defined by 
Eq. (8) differs by a factor 2 from the definition usually 
used in the literature, because of the redistribution only 
over the day side. Note also that T eq and T t herm differ 
significantly at young ages, when the intrinsic flux of the 
planet is not negligible. 

The effect of irradiation on atmosphere structures is il- 
lustrated in Fig.0for different values of the effective tem- 
perature T e ff . As already stressed in Barman et al. (2001), 
an irradiated structure characterised by Therm can differ 
significantly from a non-irradiated structure at the same 
effective temperature T e ff = Therm- This point (see also 
Seager & Sasselov 1998; Guillot & Showman 2002) empha- 
sizes the fact that adopting outer boundary conditions, 
for evolutionary calculations, from atmospheric profiles of 
nonirradiated models with T e ff = Therm, or T c ff = T cq (as 
e.g. Burrows et al. (2000)), is incorrect and yields erro- 
neous evolutionary properties for irradiated objects. 

3.2. Effect on evolution 

The main effect of irradiation on convective atmospheres 
and its consequences on evolution is well known (see 
Hubbard 1977; Brett & Smith 1993; Guillot et al 1996; 
Hubbard et al. 2002). The heating of the outer layers 
by the incident flux reduces the temperature gradient be- 
tween these layers and the interior. They become radiative 
and the top of the convective zone is displaced to larger 
depths compared to the non-irradiated clearly il- 

lustrated in Fig.El The inner atmosphere structure is hot- 
ter at a given pressure than the nonirradiated atmosphere 
model of same T e ff (see Fig- EJ - In order to match the same 
inner entropy, or the same values of P and T, character- 
istic of the boundary layer between the interior structure 
and the irradiated atmosphere structure, characterized by 
a given T e ff and log g, one would need a nonirradiated 
atmosphere model with higher T c ff, i.e a larger heat loss. 
Therefore, for a given planet heat content, i.e. internal en- 
tropy, the heat loss is reduced in the case of irradiation 
and the planet maintains a higher entropy for a longer 
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Fig. 5. Effect of irradiation on atmosphere profiles, T 
(K) versus P (dyn/cm 2 ), characterized by a surface grav- 
ity log g = 3.0, T c ff = 1000K (upper panel) and T c ff = 
100K (lower panel). Dashed lines correspond to nonir- 
radiated structures. Solid lines are irradiated models at 
a separation a = 0.046 AU from a primary with T e ff^ 
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T cq ~ 1Q30K. The squares on the curve refer to optical 
depth, defined at A = 1.2/^m, r st d = 1 and the circles to 
T st d =100. The triangles indicate the top of the convective 
zone. 



time. Since for a given mass, the interior (P,T) profile and 
thus the entropy fix the radius, the irradiated planet has 
a larger radius than the nonirradiated counterpart at a 
given time, starting from the same initial configuration. In 
other terms, gravitational contraction, which is the domi- 
nant source of energy of the planet, proceeds more slowly 
with irradiation than without it. 

Our calculations proceed as for our low-mass star or 
brown dwarf calculations, by coupling the interior and at- 
mosphere profile at a deep enough optical depth, which de- 
fines unequivocally the fundamental properties of the ob- 
ject, m, R, T c ff, L along its evolution t (Chabrier & Baraffe, 
1997). The boundary condition between inner and atmo- 
sphere structure is fixed at T stc i = 100, which corresponds 
to a range of pressure P = 0.1 - 200 bar for the whole 
range of atmosphere models 1 . The irradiated atmosphere 
models are integrated down to an optical depth T st( f = 100 



T Bt d is denned at A = 1.2/nm 
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Fig. 6. Radius (in i?j) versus age (in yr) for EGP masses 
of 1 Mj (lower panel) and 10 Mj (upper panel). The 
dashed lines are nonirradiated models. The solid lines are 
irradiated models at a separation a — 0.046 AU from the 
parent star with T e ff t = 6000K. The dotted lines are non- 
irradiated models from Burrows et al. (1997). 



for T efr > ItiOOK and r std = 10 5 for T cS < 1000K. In both 
cases, this is deep enough to reach the top of the convec- 
tive zone and to provide a good spatial resolution of these 
layers, even for the coolest models (see Fig. EJ. In any 
case, the incident flux Fi nc drops to zero at r s td << 50, 
well above the deepest layers of the atmosphere models. 
Note that for the coolest atmospheric structures, convec- 
tion does not reach the layers corresponding to r st d = 100 
(see Fig.^J. In that case, the radiative gradient in the in- 
terior is calculated with the Rosseland means of the same 
atmospheric opacities, for a consistent treatment between 
the interior and atmosphere thermal structures. 

We have calculated the evolution of planets covering 
a range of mass from 0.5 Mj to 10 Mj with and without 
irradiation. The evolution of the radius as a function of 
time is shown in Fig. for irradiated and nonirradiated 
EGPs of 1 Mj and 10 Mj. As expected, the less massive 
the planet, the larger the effect of irradiation, for a given 
incident flux. At 1 Gyr, the 0.5 Mj EGP has a 14% larger 
radius than its nonirradiated counterpart, whereas for the 
1 Mj (resp. 10 Mj), R is only 10% larger (resp. 7%). We 
also compare our COND models (the nonirradiated mod- 
els) to the Burrows et al. (1997, hereafter B97) nonirradi- 



ated models. Significant differences appear at young ages 
(< 1 Gyr), due certainly to different initial conditions (see 
Baraffe et al. 2002). For ages > 1 Gyr and m>5Mj, the 
differences between the B97 models and ours are of the 
order of the irradiation effects. This reflects the different 
input physics, mainly in the dust treatment and molec- 
ular opacities and illustrates the present uncertainties in 
the models. For ra<5Mj, however, irradiation effects be- 
come larger than the differences between the B97 and our 
models. 

The specific case of HD209458b, with a mass m p = 
0.69Mj, is illustrated in Fig. The intrinsic luminosity 
and corresponding effective temperature in the irradiated 
case (solid lines) are compared to the non-irradiated case 
(dashed lines). Starting from the same initial configuration 
in both cases, the heat loss is reduced at early ages in the 
case of irradiation, as expected. Consequently, the irradi- 
ated model evolves at larger entropy and radius than its 
non-irradiated counterpart. During the first Myr of evolu- 
tion, both evolutionary sequences contract with increasing 
central density and temperature, the non-irradiated model 
being denser. The latter becomes partially degenerate ear- 
lier, its contraction slows down and its heat loss becomes 
smaller than in the irradiated case (at logt ~ 6.2 yr). 
The situation reverses at logt ~ 7.4 yr when the effect 
of partial degeneracy becomes important in the irradiated 
sequence. The age of HD209458 is about 4-7 Gyr, accord- 
ing to Cody & Sasselov (2002). At 5 Gyr, the irradiated 
sequence displayed in Fig. predicts a radius R = 1.09i?j, 
26% smaller than the observed value. Without including 
irradiation effects, the radius is > 30% than the obser- 
vationally determined one. Note that the nonirradiated 
sequence stops at T e g = 100K, corresponding to an age of 
~ 2 Gyr and a radius R = liJj. In the following section, 
we analyse the possible reasons for such a discrepancy. 

4. Discussion 

4.1. Uncertainties of irradiated 

atmosphere/evolutionary models 

The question rises whether uncertainties of current models 
can explain the mismatch of HD209458b predicted versus 
observed radius, and whether irradiation effects can still 
provide the solution to the problem. We first note that 
our choice of parameters for the irradiated atmosphere 
calculations certainly overestimates the effects of irradia- 
tion (see §3). We assume redistribution of the incident flux 
over the day-side of the planet only, without taking into 
account varying angles of incidence of the impinging flux. 
As shown in Brett & Smith (1993), the effect of irradiation 
on the atmosphere will decrease with increasing angles of 
incidence. Note that adopting the maximum case, i.e no 
redistribution at all of -Fine, affects significantly the outer 
atmospheric profile, but only slightly the inner profile. We 
did not consider horizontal energy flows, which may ex- 
ist in a real planet receiving a varying amount of incident 
flux over its surface. But our assumption of isotropic in- 
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Fig. CJ. In order to slow down the planet contraction to 
reach the observed value at ~ 5 Gyr, the inner profiles of 
the present irradiated atmosphere models, for a given T e g, 
should be drastically modified. Such a major modification 
seems unlikely, for a fixed incident flux. Note that assum- 
ing a constant flux from the parent star during the whole 
planet evolution overestimates the total incident flux re- 
ceived by the planet over ~ 5 Gyr. For a parent star mass 
~ 1.06M© (Cody & Sasselov 2002), most of the star evo- 
lution after the first 1 Myr proceeds at a luminosity L 
smaller than its value at 5 Gyr. 

Finally, possible uncertainties due to abundance ef- 
fects, such as non-solar metallicity and/or helium abun- 
dance are not expected to affect significantly the present 
results. To estimate such uncertainties, we have computed 
a grid of irradiated atmosphere models for an over-solar 
metallicity [M/H] = +0.3 and the corresponding evolu- 
tionary sequence for a mass m p = 0.69 Mj. After a few Gyr 
evolution, the radius of the later sequence is essentially the 
same as for the solar metallicity case. Given the optimiza- 
tion of irradiation effects provided by our assumptions, 
and the huge effect required on the inner profile to reduce 
the mismatch between observed and predicted radii, we 
do not expect uncertainties in the present models to be 
the source of the discrepancy. 



Fig. 7. Effect of irradiation on the evolution of a planet 
with Trip = 0.69-Mj at a separation a = 0.046 AU from 
its parent star with T e g = 6000K. The panels from top to 
bottom display respectively the intrinsic luminosity Lint, 
the effective temperature (in K) and the radius versus time 
(in yr) . The solid curves correspond to the irradiated case 
and the dashed curves to the nonirradiated counterpart. 
We recall that in the case of irradiation, T c g and L- mt do 
not characterize the total flux emitted by the planet. The 
position of HD209458b in the lower panel is from Cody & 
Sasselov (2002). 



cident flux, with the maximum amount of flux allowed, 
should provide an upper limit for the expected effects of 
irradiation on the evolution of a planet. 

Although present atmosphere models still have some 
shortcomings, due either to incomplete molecular opaci- 
ties or to dust treatment (see §2.2), the inner profiles of 
the irradiated models must be changed and heated drasti- 
cally to provide the effect required on evolutionary models 
to match the radius of HD209458b. Fig. [7| indeed shows 
that the contraction of the irradiated sequence proceeds 
too rapidly after the first Myr of evolution to maintain a 
radius close to the observed value. At an age t\ — 1 Myr, 
the model has a radius R\ — 2Rj and an intrinsic lumi- 
nosity logLi/L = —5.12, corresponding to T e g — 670K. 
This corresponds to a characteristic Kelvin-Hclmholtz 
timescale tkh = Grn^/{R\L\) ~ lOMyr. In fact, the 
model has already reached a radius of Rio — 1.48 Rj, 
about the radius of HD209458b, after only 10 Myr (see 



4.2. Observed versus theoretical radius 

The definition of the radius in low mass stars, brown 
dwarfs or isolated giant planets is usually not a matter 
of confusion, given the negligible extension of their atmo- 
sphere compared to the total radius of the object. The 
picture could be different in the case of irradiated atmo- 
spheres, where extension effects due to the large heating of 
the upper layers may not be negligible (Seager & Sasselov 
2000; Hubbard et al. 2001). 

As discussed in Baschek et al. (1991), the condition of 
compactness for a photosphere in hydrostatic equilibrium 
is Hp/r << 1, whit Hp the pressure scale height and r 
the radial distance to the center. This condition is known 
to be perfectly fulfilled in (non-irradiated) low mass ob- 
jects (see Chabrier & Baraffe 1997), where the extension 
of the photosphere is usually less than 1% the total size 
of the object. Thus the radius can unambiguously be de- 
fined, quoting Baschek et al. (1991), as the distance of the 
atmosphere to the object center. In all our previous work, 
we fix the boundary condition between atmosphere and 
inner structure at r st d=100, knowing that i?(r st d = 100) 
is essentially the same as i?(r st d ~ 1) (see Chabrier & 
Baraffe 1997). As already mentioned, we define r s td at 1.2 
/an, which corresponds to the peak of the flux emitted by 
cool (non-irradiated) objects. Usually, the region where 
Tstd ~ 1 is close to the region where rR OSSe iand ~ 1 (for the 
present irradiated models as well). Above this region, the 
atmosphere contains little mass and contributes negligibly 
to the luminosity. The evolutionary calculations for irradi- 
ated models presented in §3.2 determine also the radius at 
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T st d = 100. Within the present assumptions of irradiation, 
the atmospheric extension between T st d=100 and r s td ~ 1 
represents only 1-2% of the total radius. The presently 
calculated theoretical radius is thus essentially equivalent 
to a photospheric radius at 1.2 fim. 

However, the observed radius of HD209458b (e.g the 
one estimated by Cody & Sasselov 2002) is based on the 
analysis of optical light curves. It corresponds to a region 
of the atmosphere where the optical depth is near unity 
at 0.6/xm, which is near r st d = 10~ 2 in our atmosphere 
models. Therefore, the radius predicted by the evolution- 
ary models is not equivalent to the observed radius. If the 
atmospheric extension and the opacity of the atmosphere 
are large at the observed wavelengths, the measured ra- 
dius could be different from the radius predicted by the 
evolutionary calculation. However, based on our irradiated 
atmosphere models with the gravity predicted by the evo- 
lution and the age of HD209458b (i.e. log 5=3.2), the at- 
mospheric extension between r st d = 100 and r s td = 10 -2 
(where the optical depth is close to unity at 0.6/1 m) is 
very small, namely 0.05 compared to the overall ra- 
dius (Barman et al. 2003, in preparation). Adding this 
value to the radius predicted by the models at 5 Gyr 
yields an optical-depth radius at 0.6/xm of ~ 1.14 still 
22% less than the observed value. For younger planets or 
planets undergoing stronger irradiation effects, with much 
lower gravities, the extension will be more important and 
should be taken into account for a consistent comparison 
between theoretical and observed radii, as already stressed 
by Seager & Sasselov (2000) and Hubbard et al. (2001). 

In the same vein, Hubbard et al. (2001) estimate a 
radius of 94 430 km (1.32 i?j) at a pressure of 1 bar, based 
on a detailed analysis of physical effects influencing the 
observed light curve of HD 209458b. At 5 Gyr, our models 
predict a radius at 1 bar of ~ 1.1 Rj, 18% smaller than 
the Hubbard et al. (2001) estimate. Such a discrepancy 
is consistent with the afore-mentioned mismatch for the 
radius at 0.6/xm. 

4.3. Other sources of energy deposit 

If irradiation effects alone do not explain the large ob- 
served radius of HD 209458b, other sources of energy 
must be invoked. Tidal interactions between the star 
and the planet can provide a source of energy associ- 
ated to the synchronization and/or circularization of the 
planet orbit, dissipated within the planet (Lubow et al. 
1997; Rieutord & Zahn 1997; Bodenheimer et al. 2001). 
However, as discussed recently by Guillot & Showman 
(2002) and Showman & Guillot (2002), these processes 
are efficient only during the early stages of the planet 
evolution. Estimates based on the current understand- 
ing of such processes yield typical circularisation timescale 
Tcirc ~ 10 8 yr (Bodenheimer et al. 2001) and synchronisa- 
tion timescale r syn < 10 s yr (Lubow et al. 1997; Rieutord 
& Zahn 1997). Such an energy source seems thus unlikely 
to slow down the long term evolution of the planet, un- 



less a second planet orbiting HD209458a is present. Such 
a detection has been claimed very recently in the litera- 
ture (Bodenheimer, Laughlin & Lin, 2003) but remains to 
be confirmed unambiguously. Showman & Guillot (2002) 
suggested that downward transport of kinetic energy pro- 
duced by atmospheric circulation could be dissipated in 
the planet interior, leading to a substantial deposit of 
energy. Within the present input physics and treatment 
of irradiation, we can estimate the amount of energy re- 
quired to reach the radius of HD 209458b. As in Guillot 
& Showman (2002), we arbitrarily add an extra term of 
energy generation e cx tra in the energy equation at different 
depths. We have explored several cases displayed in Fig.|SJ 
We add a total amount L G xtra = / ^cxtradm in a region of 
mass Am enclosed between the surface and an arbitrary 
depth at mass shell mi (i.e Am = m p —mi). Various tests 
indicate that an amount of energy L oxtra ~ 10 27 — 5 x 10 27 
erg.s -1 dissipated along the internal adiabat yields a ra- 
dius within the error bars of the observed value (see Fig. 
IHJ) . As expected, the larger the fraction of L e xtra deposited 
in the convective layers, the more important the effect. 
Note that the case displayed in Fig. with L e xtra ~ 10 27 
erg.s" 1 dissipated all over the star (Am/m p — 1, dashed 
curve) is equivalent to depositing the same amount of en- 
ergy only at the very center. Our quantitative estimates 
are in general agreement with Guillot & Showman (2002). 
Such an amount of energy represents more than 100 times 
the intrinsic luminosity L; nf of the planet, which is ~ 10 25 
erg.s -1 at 1 Gyr and ~ 2 1 24 erg.s^ 1 at 5 Gyr (see Fig. 
EJ. However, it represents only ~ 1% of the incoming lu- 
minosity, L inc = 2nRpF inc ~ 10 29 erg s _1 , which largely 
dominates the planet energetic balance. Thus, an alterna- 
tive possibility is the release of an external source of en- 
ergy caused by the incident radiation (see e.g. Showman 
& Guillot, 2002). As illustrated in Fig. 03 however, the 
extra source of energy must be dissipated at the top of 
the internal adiabat, i.e. at a much deeper level than the 
penetration of the incident photons (>> r = 1). 

5. Conclusion 

We have presented calculations describing the evolution 
of cool brown dwarfs and extra-solar gaseous planets. The 
present models reproduce the main trends of observed 
methane-dwarfs in near-IR color- magnitude diagrams ( J- 
K, K-L). Problems still remain at wavelengths < 1 /xm, 
with a flux excess predicted in the I-bandpass. The treat- 
ment of atomic line broadening in such dense objects may 
be the source of the present discrepancy. The models fail 
to reproduce the coolest L-dwarfs and a detailed treat- 
ment of dust diffusion/sedimentation is required for a cor- 
rect description of the transition region between L- and 
T- dwarfs. Work is in progress in this direction. 

We have included the effects of irradiation, coupling ir- 
radiated atmosphere profile and inner structure, and pro- 
viding consistent evolutionary models for irradiated plan- 
ets. The effect of irradiation are shown to modify signif- 
icantly the mechanical (mass-radius) and thermal evolu- 
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Fig. 8. Effect of extra source of energy dissipation on 
the evolution of a planet with mass m p = 0.69Mj. The 
solid and long-dashed lines correspond to the irradiated 
and non-irradiated sequences respectively, with no extra 
source of energy dissipation. The other curves correspond 
to irradiated sequences with a total amount of additional 
energy L e xtra (i n erg.s" 1 ) deposited in layers between the 
surface and the mass shell m p — Am, as indicated on the 
figure (see text). 



tion of irradiated EGPs. However, a significant discrep- 
ancy (26%) remains between the theoretical and observed 
radii of the transit planet HD209458b. We have explored 
possible uncertainties inherent in the models to explain 
such a discrepancy. Although solving these uncertainties 
may modify the outer structure of the models (extension 
of the atmosphere, albedo) and perhaps slightly reduce the 
discrepancy, none of the uncertainties is likely to modify 
significantly the inner entropy profile of the models, which 
determines the radius of the planet. Indeed, a drastic mod- 
ification on the inner thermal structure is required to bring 
the theoretical radius in agreement with the observed one. 

In summary, we do not expect irradiation effects alone 
to explain the large observed radius of HD209458b. In 
the same vein, tidal interactions will affect only the early 
stages of evolution of the planet but will probably be dis- 
sipated too rapidly to affect the long term contraction of 
the object. Other sources of energy, representing about 
100 times the intrinsic luminosity of the planet, seem 
to be required to explain the observed radius. The first 



extra-solar planet transit thus remains a challenge for the- 
ory. Detection of other transits is now crucial to conclude 
whether HD209458 is a peculiar system, whether a sec- 
ond planetary companion is confirmed or not, or whether 
we are missing something in the current understanding of 
close-in giant planets. 

Note: Isochrones for t > 1 Myr of the COND models 
(from 0.5 Mj to 0.1 Mq) are available at: 
http://www.ens-lyon.fr/" ibaraffe/COND03_models 
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Table 1. COND isochroncs for 0.1 Gyr 
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16.16 


16.01 


16.38 


13.26 


12.55 





.0080 


816. 


-5. 


.246 


0. 


120 


4. 


180 


26 


.77 


22 


,28 


19. 


,49 


15.76 


15.65 


15.79 


12.97 


12.35 





.0090 


886. 


-5. 


.103 


0. 


120 


4, 


232 


26 


.45 


21 


,96 


19 


,19 


15.32 


15.23 


15.16 


12.63 


12.13 





.0100 


953. 


-4 


.978 


0. 


120 


4. 


279 


26 


.10 


21. 


,66 


18, 


,92 


14.94 


14.86 


14.69 


12.34 


11.96 





.0120 


1335. 


-4 


.332 


0. 


129 


4. 


.297 


23 


.53 


19 


,44 


16, 


,79 


13.20 


12.97 


12.76 


10.90 


11.17 





.0150 


1399. 


-4 


.281 


0. 


124 


4. 


424 


23 


.30 


19. 


,24 


16 


.46 


13.05 


12.82 


12.65 


10.83 


11.15 





.0200 


1561. 


-4. 


.110 


0. 


122 


4, 


569 


22 


.30 


18. 


,55 


16. 


,08 


12.60 


12.34 


12.17 


10.53 


10.99 





.0300 


1979. 


-3 


.668 


0. 


126 


4, 


715 


19 


.96 


16. 


,80 


14. 


.48 


11.52 


11.20 


10.90 


9.82 


10.38 





.0400 


2270. 


-3 


.386 


0. 


132 


4. 


797 


18. 


.46 


15. 


.63 


13 


.31 


10.89 


10.52 


10.19 


9.39 


9.84 





.0500 


2493. 


-3 


.167 


0. 


141 


4. 


837 


17 


.09 


14. 


.77 


12 


,53 


10.43 


10.02 


9.71 


9.04 


9.37 





.0600 


2648. 


-3 


.008 


0. 


150 


4. 


863 


16 


.08 


14. 


,12 


12 


.01 


10.10 


9.68 


9.37 


8.78 


9.03 





.0700 


2762. 


-2 


.879 


0. 


160 


4. 


874 


15 


.33 


13 


,59 


11 


.60 


9.82 


9.39 


9.10 


8.55 


8.75 





.0720 


2782. 


-2 


.856 


0. 


162 


4, 


875 


15 


.20 


13. 


,50 


11 


.53 


9.77 


9.34 


9.05 


8.51 


8.70 





.0750 


2809. 


-2 


.821 


0. 


166 


4. 


875 


15 


.01 


13. 


36 


11 


.42 


9.69 


9.26 


8.97 


8.44 


8.63 





.0800 


2846. 


-2 


.776 


0. 


170 


4. 


880 


14 


.77 


13. 


,18 


11. 


.29 


9.60 


9.16 


8.87 


8.36 


8.53 





.0900 


2910. 


-2 


.689 


0. 


180 


4. 


884 


14 


.34 


12. 


,85 


11 


.03 


9.40 


8.96 


8.68 


8.19 


8.35 





.1000 


2960. 


-2 


.617 


0. 


189 


4. 


887 


14 


.02 


12 


,58 


10, 


,82 


9.24 


8.80 


8.52 


8.05 


8.19 



Notes.- T e ff is in K, the gravity g in cgs. The VRI magnitudes are in the Johnson- Cousins system (Bessell 1990), JHK in 
the CIT system (Leggett 1992), L' in the Johnson-Glass system and M in the Johnson system. 



Table 2. Same as Table 1 for 0.5 Gyr 



m/M© 


T eff 


log L/L e 


R/Rq 


log g 


M v 


Mr 


Mi 


Mj 


M H 


M K 


M v 


M M 


0.0005 


141. 


-8.415 


0. 


105 


3.097 


56.30 


51.03 


46.60 


37.42 


33.07 


51.62 


23.09 


20.59 


0.0010 


203. 


-7.753 


0. 


109 


3.365 


47.57 


42.88 


38.99 


31.61 


29.15 


43.23 


20.93 


18.68 


0.0020 


272. 


-7.218 


0. 


112 


3.639 


37.05 


33.00 


30.06 


25.07 


24.62 


34.02 


18.66 


16.58 


0.0030 


322. 


-6.913 


0. 


113 


3.805 


32.02 


28.23 


25.75 


22.05 


22.27 


29.03 


17.52 


15.53 


0.0040 


370. 


-6.670 


0. 


114 


3.928 


30.65 


26.73 


24.16 


21.01 


21.06 


26.45 


16.85 


15.05 


0.0050 


409. 


-6.496 


0. 


113 


4.027 


29.60 


25.57 


22.94 


20.20 


20.11 


24.54 


16.32 


14.67 


0.0060 


449. 


-6.340 


0. 


113 


4.112 


29.16 


25.05 


22.39 


19.64 


19.51 


23.10 


15.92 


14.36 


0.0070 


488. 


-6.200 


0. 


112 


4.185 


28.71 


24.51 


21.80 


19.10 


18.91 


21.68 


15.52 


14.06 


0.0080 


525. 


-6.080 


0. 


111 


4.249 


28.40 


24.14 


21.41 


18.65 


18.46 


20.74 


15.19 


13.83 


0.0090 


564. 


-5.963 


0. 


110 


4.307 


28.14 


23.82 


21.07 


18.21 


18.04 


19.95 


14.88 


13.63 


0.0100 


599. 


-5.864 


0. 


110 


4.358 


27.91 


23.53 


20.77 


17.80 


17.66 


19.22 


14.59 


13.43 


0.0120 


759. 


-5.447 


0. 


110 


4.432 


27.20 


22.73 


19.91 


16.38 


16.29 


16.82 


13.51 


12.68 


0.0150 


791. 


-5.404 


0. 


107 


4.557 


27.11 


22.63 


19.81 


16.20 


16.16 


16.58 


13.41 


12.60 


0.0200 


936. 


-5.133 


0. 


104 


4.704 


26.53 


22.07 


19.27 


15.33 


15.34 


15.37 


12.78 


12.21 


0.0300 


1264. 


-4.636 


0. 


101 


4.905 


24.97 


20.77 


17.95 


13.90 


13.87 


13.69 


11.70 


11.68 


0.0400 


1583. 


-4.255 


0. 


100 


5.040 


23.11 


19.28 


16.74 


12.92 


12.76 


12.68 


10.95 


11.33 


0.0500 


1875. 


-3.955 


0. 


101 


5.131 


21.31 


17.94 


15.53 


12.21 


11.95 


11.79 


10.44 


10.99 


0.0600 


2116. 


-3.729 


0. 


102 


5.194 


19.99 


16.91 


14.50 


11.69 


11.36 


11.13 


10.10 


10.63 


0.0700 


2329. 


-3.534 


0. 


106 


5.233 


18.84 


16.04 


13.68 


11.27 


10.90 


10.63 


9.81 


10.26 


0.0720 


2369. 


-3.498 


0. 


107 


5.238 


18.60 


15.89 


13.54 


11.19 


10.81 


10.54 


9.75 


10.19 


0.0750 


2426. 


-3.445 


0. 


108 


5.244 


18.25 


15.66 


13.33 


11.08 


10.69 


10.42 


9.67 


10.08 


0.0800 


2518. 


-3.356 


0. 


111 


5.248 


17.65 


15.27 


13.00 


10.89 


10.49 


10.22 


9.53 


9.89 


0.0900 


2680. 


-3.189 


0. 


119 


5.241 


16.54 


14.56 


12.43 


10.55 


10.12 


9.85 


9.25 


9.53 


0.1000 


2804. 


-3.047 


0. 


128 


5.223 


15.68 


13.96 


11.98 


10.25 


9.80 


9.54 


9.00 


9.22 



Baraffe et al.: Evolutionary models for cool brown dwarfs 
Table 3. Same as Table 1 for 1 Gyr 



ra/M© T eff log L/Lq R/Rq log g M v M R Mi Mj M H M K M L , M M 






.0005 


111. 


-8 


.851 


0. 


102 


3 


.115 


60. 


75 


55 


23 


50 


.50 


40 


.19 


35.07 


55 


.87 


24.15 


21.49 





.0010 


160. 


-8. 


.185 


0. 


106 


3 


.386 


54. 


15 


49. 


.10 


44, 


.69 


35 


.58 


32.06 


49 


.17 


22.40 


19.95 





.0020 


226. 


-7 


.560 


0. 


109 


3 


.662 


44. 


39 


39 


,91 


36 


.34 


29 


.28 


27.80 


40. 


.31 


20.18 


17.94 





.0030 


270. 


-7. 


.244 


0. 


111 


3 


.827 


37, 


64 


33 


,60 


30. 


.73 


25 


.29 


24.99 


34. 


.51 


18.84 


16.66 





.0040 


304. 


-7. 


.031 


0. 


111 


3 


.950 


32, 


62 


28. 


,93 


26. 


.60 


22 


.49 


22.91 


30 


.25 


17.91 


15.73 





.0050 


342. 


-6 


.831 


0. 


110 


4 


.051 


31, 


58 


27 


,79 


25. 


.36 


21 


.71 


21.98 


28. 


.28 


17.38 


15.36 





.0060 


377. 


-6 


.664 


0. 


110 


4 


.134 


30, 


53 


26 


63 


24 


.12 


20 


.96 


21.07 


26. 


.39 


16.87 


15.01 





.0070 


403. 


-6 


.556 


0. 


109 


4 


.208 


29, 


.77 


25. 


,79 


23 


.26 


20 


.41 


20.43 


25. 


.04 


16.54 


14.76 





.0080 


438. 


-6 


.417 


0. 


108 


4 


.272 


29, 


37 


25. 


,31 


22. 


.73 


19 


.93 


19.89 


23. 


.77 


16.18 


14.50 





.0090 


464. 


-6 


.325 


0. 


107 


4 


.331 


29, 


.06 


24 


,94 


22. 


.33 


19 


.59 


19.49 


22. 


.85 


15.93 


14.31 





.0100 


491. 


-6 


.235 


0. 


107 


4 


.383 


28, 


,74 


24. 


.55 


21. 


.89 


19 


.23 


19.07 


21. 


.90 


15.65 


14.11 





.0120 


578. 


-5 


.955 


0. 


106 


4 


.467 


28, 


,09 


23 


,75 


21. 


.01 


18. 


.15 


18.03 


19 


.86 


14.88 


13.60 





.0150 


628. 


-5 


.835 


0. 


103 


4 


.587 


27. 


86 


23 


,47 


20 


.72 


17 


.70 


17.62 


19. 


.05 


14.56 


13.38 





.0200 


766. 


-5. 


.514 


0. 


100 


4 


.736 


27, 


31 


22. 


,85 


20. 


.05 


16 


.56 


16.55 


17 


.19 


13.72 


12.80 





.0300 


1009. 


-5. 


.071 


0. 


096 


4. 


.948 


26, 


.40 


21. 


,96 


19. 


.15 


15 


.10 


15.16 


15 


.14 


12.67 


12.15 





.0400 


1271. 


-4 


.696 


0. 


093 


5 


.099 


25, 


19 


20 


,99 


18. 


.13 


14 


,04 


14.04 


13 


.90 


11.88 


11.80 





.0500 


1543. 


-4 


.374 


0. 


092 


5 


.211 


23, 


73 


19 


,81 


17 


.15 


13 


.21 


13.12 


13 


.04 


11.25 


11.53 





.0600 


1801. 


-4 


.106 


0. 


092 


5 


.291 


22, 


,13 


18 


,59 


16 


.10 


12 


.56 


12.36 


12 


.27 


10.77 


11.26 





.0700 


2082. 


-3 


.829 


0. 


094 


5 


.333 


20, 


,44 


17 


.31 


14 


.87 


11 


.93 


11.62 


11 


.43 


10.32 


10.85 





.0720 


2140. 


-3 


.772 


0. 


095 


5 


.336 


20, 


.12 


17. 


.05 


14 


.61 


11 


.80 


11.48 


11 


.27 


10.23 


10.75 





.0750 


2234. 


-3. 


.679 


0. 


098 


5 


.334 


19, 


59 


16, 


.63 


14. 


.21 


11 


.60 


11.25 


11 


.02 


10.08 


10.58 





.0800 


2383. 


-3 


.527 


0. 


102 


5 


.323 


18, 


.67 


15, 


.96 


13. 


.61 


11 


.26 


10.89 


10 


.63 


9.84 


10.27 





.0900 


2627. 


-3. 


.268 


0. 


113 


5 


.285 


16. 


,98 


14 


,86 


12 


.67 


10 


.72 


10.30 


10 


.03 


9.40 


9.72 





.1000 


2784. 


-3. 


.083 


0. 


125 


5 


.246 


15, 


86 


14. 


,09 


12 


.09 


10 


.33 


9.89 


9.62 


9.07 


9.31 



Table 4. Same as Table 1 for 5 Gyr 



m/M Q T eff log L/Lq R/Rq log g M v M R Mj Mj M H M K M v M M 






.0020 


129. 


-8. 


.570 


0. 


.105 


3 


.698 


60 


.05 


54. 


.63 


49. 


.68 


38 


.16 


34.52 


53.62 


23.33 


20.80 





.0030 


162. 


-8. 


.166 





.105 


3 


.868 


55 


.32 


50. 


.15 


45. 


.63 


34 


.96 


32.46 


49.21 


22.22 


19.81 





.0040 


193. 


-7 


.867 





.105 


3 


.994 


50 


.68 


45 


,80 


41 


.76 


32 


.14 


30.52 


45.15 


21.28 


18.92 





.0050 


220. 


-7, 


.644 





.105 


4 


.095 


46 


.50 


41 


,92 


38. 


.30 


29 


.83 


28.84 


41.71 


20.48 


18.16 





.0060 


244. 


-7. 


.469 





.104 


4 


.180 


42 


.71 


38. 


,43 


35. 


.18 


27 


.80 


27.34 


38.67 


19.78 


17.49 





.0070 


265. 


-7. 


.328 





.103 


4 


.254 


39 


.29 


35 


,29 


32. 


.37 


26 


.02 


26.00 


35.96 


19.16 


16.88 





.0080 


284. 


-7. 


.217 





.103 


4 


.318 


36 


.31 


32. 


.57 


29. 


.93 


24 


.50 


24.85 


33.65 


18.64 


16.36 





.0090 


301. 


-7 


.124 





.102 


4 


.376 


33 


.73 


30 


,22 


27 


.82 


23 


.21 


23.85 


31.64 


18.18 


15.91 





.0100 


322. 


-7 


.015 





.101 


4 


.429 


33 


.05 


29 


,47 


27 


.06 


22 


.74 


23.30 


30.49 


17.91 


15.71 





.0120 


361. 


-6 


.823 





.100 


4 


.519 


31 


.58 


27 


,86 


25 


.46 


21 


.77 


22.17 


28.16 


17.41 


15.32 





.0150 


399. 


-6 


.671 





.098 


4 


.634 


30 


.20 


26. 


,35 


24. 


.00 


20 


.85 


21.14 


26.01 


16.96 


14.97 





.0200 


473. 


-6. 


.401 





.095 


4 


.786 


29 


.24 


25. 


,17 


22 


.66 


19 


.89 


19.92 


23.31 


16.18 


14.45 





.0300 


610. 


-6. 


.008 





.090 


5 


.011 


28 


.17 


23. 


,83 


21. 


.17 


18 


.34 


18.31 


20.04 


15.04 


13.71 





.0400 


760. 


-5. 


.670 





.085 


5 


.179 


27 


.58 


23. 


,15 


20. 


.40 


17 


.04 


17.11 


17.99 


14.17 


13.13 





.0500 


931. 


-5 


.353 





.082 


5 


.313 


27. 


.09 


22. 


.63 


19. 


.82 


15 


.94 


16.05 


16.38 


13.38 


12.63 





.0600 


1120. 


-5. 


.058 





.079 


5 


.418 


26 


.44 


22 


,03 


19 


.20 


15 


.01 


15.13 


15.15 


12.73 


12.27 





.0700 


1524. 


-4 


.504 





.081 


5 


.466 


24 


.33 


20. 


,36 


17. 


.60 


13 


.52 


13.50 


13.44 


11.60 


11.77 





.0720 


1712. 


-4 


.278 





.083 


5 


.453 


23 


.11 


19. 


,36 


16. 


.75 


12 


.97 


12.85 


12.80 


11.16 


11.54 





.0750 


2006. 


-3 


.942 





.089 


5 


.411 


21 


.03 


17. 


,80 


15. 


.34 


12 


.20 


11.92 


11.78 


10.54 


11.06 





.0800 


2320. 


-3 


.603 





.099 


5 


.353 


19 


.11 


16. 


,28 


13 


.89 


11 


.44 


11.07 


10.82 


9.97 


10.43 





.0900 


2622. 


-3 


.275 





.113 


5 


.289 


17 


.02 


14. 


,88 


12 


.70 


10 


.73 


10.31 


10.04 


9.42 


9.73 





.1000 


2785. 


-3 


.083 





.125 


5 


.247 


15 


.85 


14 


,09 


12 


.09 


10 


.33 


9.89 


9.62 


9.07 


9.31 
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Table 5. Same as Table 1 for 10 Gyr 



m/M© 


T eff 


log L/Lq 


R/Rq 


log g 


M v 


Mr 


Mi 


Mj 


M H 


M K 


M L , 


M M 


n nnon 
U.UUoU 


125. 


q c?on 
-o.Dza 


U.1U4 


O.S (9 


c? i /if? 

D1.4D 


on 

55.89 


CA 7/ A 

5u. / y 


OQ OA 
OO.Z9 


O C A/1 

OO.U4 


C /I OQ 

54.28 


OO OC? 


oa on 
2U.89 


n fin a n 
U.UU4U 


149. 


o ooc 
-8.325 


a i nA 
U.104 


/I A AC? 

4.UUb 


C O AA 

58. UU 


co c?o 
52. Dz 


<7 oc 

4/.8b 


OC? Al 

36.01 


OO cc? 

33.56 


CI 11 

51.11 


OO C A 

22.59 


OA 1 A 

2U.19 


U.UUoU 


1 IZ. 


-8.08/ 


A 1 AO 

U.103 


A 1 AA 

4.1U9 


c a c?o 
54. DZ 


A A /I C? 

49.4b 


A C AC 

45.05 


O A AC? 

34. U6 


OO 1 o 

32.18 


/I Q OC? 

48.2b 


1 AO 

21.92 


1 A C T 

19.5 1 


n nni^n 
U.UUoU 


i no 
193. 


T OOO 

- ( .888 


A 1 AO 

0.102 


/I 1 AC 

4.195 


CI OA 

51.29 


A C? OO 

4b. 38 


/I O OA 

42.30 


OO OO 

61.16 


OA Q C 

30.85 


A C CO 

45.53 


1 OA 

21.29 


1 O AT 

18.9 1 


U.UU fU 


Ol o 

213. 


T 70 /I 

-f ./24 


A 1 AO 

0.102 


/I OTA 

4.2 ( U 


AO IT 

48.1 ( 


/i o 

43.50 


OA TO 

39. (1 


OA C C 

30.55 


OA C?0 

2y.o2 


/I O Al 

43.01 


OA T1 

2U. 1 1 


1 O A 1 

18.41 


U.UUoU 


OQO 

232. 


T C O A 

- 1 .584 


A 1 Al 
U.101 


A OOC 

4.335 


45.19 


/I A 77 

4U. / / 


OT OO 

61 .28 


OA AA 

29. uo 


OO AC 

28.4b 


/I A C?C? 

4U.bb 


OA 1 Q 

2U. 18 


1 T OA 

1 / .89 


u.uuyu 


249. 


t /i c? a 


A 1 AA 
U.100 


A 9 AO 

4.393 


/I O A A 

42.49 


OO OA 

38.30 
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